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Abstract 

We compute metric correlations in loop quantum gravity with the dynamics defined by the 
new spin foam models. The analysis is done at the lowest order in a vertex expansion and at 
the leading order in a large spin expansion. The result is compared to the graviton propagator 
of perturbative quantum gravity. 

1 Introduction 

In this paper we compute metric correlations in Loop Quantum Gravity (LQG) [HOE] and we 
compare them with the scaling and the tensorial structure of the graviton propagator in perturbative 
Quantum Gravity [H El E] . The strategy is the one introduced in [7] and developed in [U EJ [101 
[TT] [T2l [T3| Q31 02]. In particular, we use the boundary amplitude formalism [TJ [16j [T71 [18]. The 
dynamics is implemented in terms of (the group field theory expansion of) the new spin foam 
models introduced by Engle, Pereira, Rovelli and Livine (EPRL 7 model) [19] and by Ereidel and 
Krasnov (FK 7 model) [20] . We restrict attention to Euclidean signature and Immirzi parameter 
smaller than one: < 7 < 1. In this case the two models coincide. 

Previous attempts to derive the graviton propagator from LQG adopted the Barrett-Crane spin 
foam vertex (21] as model for the dynamics [3 El El [IH1 El 121 H3 El [IS] (see also [22l [23l [24] 
for investigations in the three-dimensional case) . The analysis of |121 [T3] shows that the Barrett- 
Crane model fails to give the correct scaling behavior for off-diagonal components of the graviton 
propagator. The problem can be traced back to a missing coherent cancellation of phases between 
the intertwiner wave function of the semiclassical boundary state and the intertwiner dependence of 
the model. The attempt to correct this problem was part of the motivation for the lively search of 
new spin foam models with non-trivial intertwiner dependence |25 [ \26 \ \19 \ \20 \ I27j. The intertwiner 
dynamics of the new models was investigated numerically in [281 l29"l I5UI I5T] . The analysis of the 
large spin asymptotics of the vertex amplitude of the new models was performed in 
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and in [35]. In [14j . the obstacle that prevented the Barrett-Crane model from yielding the correct 
behaviour of the propagator was shown to be absent for the new models: the new spin foams 
feature the correct dependence on intertwiners to allow a coherent cancellation of phases with the 
boundary semiclassical state. In this paper we restart from scratch the calculation and derive the 
graviton propagator from the new spin foam models. 

In this introduction we briefly describe the quantity we want to compute. We consider a manifold 
1Z with the topology of a 4-ball. Its boundary is a 3-manifold £ with the topology of a 3-sphere 
S 3 . We associate to £ a boundary Hilbert space of states: the LQG Hilbert space TLy, spanned by 
(abstract) spin networks. We call l^) a generic state in TIy,- A spin foam model for the region 1Z 
provides a map from the boundary Hilbert space to C. We call this map {W\. It provides a sum 
over the bulk geometries with a weight that defines our model for quantum gravity. The dynamical 
expectation value of an operator O on the state |\&) is defined via the following expression 1 

{0) (W\0\9) 



(W|tf) 



The operator O can be a geometric operator as the area, the volume or the length [361 EH EH 
[3U], |4"01 |4"T| I42j . The geometric operator we are interested in here is the (density-two inverse-) 
metric operator q ab (x) = 5 l i Ef(x)Ej(x). We focus on the connected two-point correlation function 
G abcd {x,y) on a semiclassical boundary state |^o)- It is defined as 

G abcd (x, y) = (q ab {x) q cd {y)) - (q ab (x)) (q cd (y)) . (3) 

The boundary state |\Po) is semiclassical in the following sense: it is peaked on a given configuration 
of the intrinsic and the extrinsic geometry of the boundary manifold S. In terms of Ashtekar- 
Barbero variables these boundary data correspond to a couple (Eq,Aq). The boundary data are 
chosen so that there is a solution of Einstein equations in the bulk which induces them on the 
boundary. A spin foam model has good semiclassical properties if the dominant contribution to the 
amplitude (ly^o) comes from the bulk configurations close to the classical 4-geometries compatible 
with the boundary data (Eq,Ao). By classical we mean that they satisfy Einstein equations. 

The classical bulk configuration we focus on is flat space. The boundary configuration that we 
consider is the following: we decompose the boundary manifold S 3 in five tetrahedral regions with 
the same connectivity as the boundary of a 4-simplex; then we choose the intrinsic and the extrinsic 
geometry to be the ones proper of the boundary of a Euclidean 4-simplex. By construction, these 
boundary data are compatible with flat space being a classical solution in the bulk. 

For our choice of boundary configuration, the dominant contribution to the amplitude (Wl^o) is 
required to come from bulk configurations close to flat space. The connected two-point correlation 



1 This expression corresponds to the standard definition in (perturbative) quantum field theory where the vacuum 
expectation value of a product of local observables is defined as 

/ D[<p]0{xi) ■ ■ ■ 0{x n )e slv] [ D[4>]W[4>]0{xi) ■ ■ ■ O(x n )^ [4>] 
(O(xi)---O(x n ))o = * - r = J - ? • (1) 



The vacuum state 9o[<f>] codes the boundary conditions at infinity. 
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function G (x, y) probes the fluctuations of the geometry around the classical configuration given 
by flat space. As a result it can be compared to the graviton propagator computed in perturbative 
quantum gravity. 

The plan of the paper is the following: in section [2] we introduce the metric operator and 
construct a semiclassical boundary state; in section [3] we recall the form of the new spin foam 
models; in section H] we define the LQG propagator and provide an integral formula for it at the 
lowest order in a vertex expansion; in section [5] we compute its large spin asymptotics; in section 
[6] we discuss expectation values of metric operators; in section [7J we present our main result: the 
scaling and the tensorial structure of the LQG propagator at the leading order of our expansion; in 
section [8] we attempt a comparison with the graviton propagator of perturbative quantum gravity. 

2 Semiclassical boundary state and the metric operator 

Semiclassical boundary states are a key ingredient in the definition of boundary amplitudes. Here 
we describe in detail the construction of a boundary state peaked on the intrinsic and the extrinsic 
geometry of the boundary of a Euclidean 4-simplex. The construction is new: it uses the coherent 
intertwiners of Livine and Speziale [27] (see also [M]) together with a superposition over spins as 
done in [8] . It can be considered as an improvement of the boundary state used in [121 El EI] 
where Rovelli-Speziale gaussian states [13] for intertwiners were used. 

We consider a simplicial decomposition A5 of S 3 . The decomposition A5 is homeomorphic to 
the boundary of a 4-simplex: it consists of five cells t a which meet at ten faces f a b (a,b = 1, . . , 5 
and a < b). Then we consider the sector of the Hilbert space spanned by spin network states 
with graph T5 dual to the decomposition A5, 



v 5 I45 V4 




T5 is a complete graph with five nodes. We call v a its nodes and l a b (a < b) its ten links. Spin 
network states supported on this graph are labelled by ten spins jab ( a < b) and five intertwiners 
i a . We denote them by \T5,j a b,i a ) and call Tir 5 the Hilbert space they span. On Hr 5 we can 
introduce a metric operator smearing the electric field on surfaces dual to links, i.e. considering 
scalar products of fluxes. We focus on the node n and consider a surface f na which cuts the link 
from the node n to the node a. The flux operator through the surface f na , parallel transported in 
the node n, is denoted 2 (-E^)i- It has the following three non-trivial properties: 
(i) the flux operators (E£)i and are related by a SU(2) parallel transport g an from the node 

2 Throughout the paper i, j, k ... = 1,2,3 are indices for vectors in R 3 . 
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a to the node n together with a change of sign which takes into account the different orientation 
of the face f an , 

(K)i = -(Ran)i j (E'a)j , (5) 

where R an is the rotation which corresponds to the group element g an associated to the link l an , 
i.e. Ran = DW{g ab )- 

(ii) the commutator of two flux operators for the same face f na is 3 

[{E a n ) i ,{E a n ) j ]=i 1 e ij k {K)k\ (6) 

(iii) a spin network state is annihilated by the sum of the flux operators over the faces bounding a 
node 

J2(KW5jah,ia)=0. (7) 

This last property follows from the SU(2) gauge invariance of the spin network node. 

Using the flux operator we can introduce the density-two inverse-metric operator at the node n, 
projected in the directions normal to the faces f na and j nb . It is defined as E®-E^ = 5^ (E^)i(E^)j. 
Its diagonal components E^-E^ measure the area square of the face f na , 

E%-EZ\T 5 ,j ab ,i a ) = (7Vjna(jna + l)) 2 \T 5 ,jab,ia) ■ (8) 

Spin network states are eigenstates of the diagonal components of the metric operator. On the 
other hand, the off-diagonal components E^-E b n with a ^ b measure the dihedral angle between 
the faces f na and /„;, (weighted with their areas). It reproduces the angle operator [40]. Using the 
recoupling basis for intertwiner space, we have that in general the off-diagonal components of the 
metric operator have non-trivial matrix elements 

K' E n |r 5 ,jaMa) = E ( E n~ E n) ^ FsJabJa) ■ (9) 

We refer to [121 113] for a detailed discussion. In particular, from property (ii), we have that some 
off-diagonal components of the metric operator at a node do not commute [39] 

From this non-commutativity an Heisenberg inequality for dispersions of metric operators follows. 
Here we are interested in states which are peaked on a given value of all the off-diagonal components 
of the metric operator and which have dispersion of the order of Heisenberg's bound. Such states 
can be introduced using the technique of coherent intertwiners [27^ 134] . A coherent intertwiner 
between the representations ji, . . , j& is defined as 4 



1 r 4 

^-^(n 1 ,..,n 4 )= = / dh Y[(ja,m a \D^\h)\j a ,n a ) (11 

VS2(m, . . ,n 4 ) Jsu(2) ^ 



3 Throughout the paper we put c = h = GNcwton = 1. 

4 The state \j, ft) is a spin coherent state. It is labelled by a unit vector ft or equivalently by a point on the unit 
sphere. Given a SI/ (2) transformation g which acts on the vector +e z sending it to the vector n = Re z , a spin 
coherent state is given by \j, n) = D^'(g)\j,+j). As a result, it is defined up to a phase e %a ^ corresponding to a 
transformation exp(ian ■ J) \j, n) — e la ^ \j, ft) . A phase ambiguity in the definition of the coherent intertwiner (|lll) 
follows. Such ambiguity becomes observable when a superposition over j is considered. 
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and is labelled by four unit vectors Hi, . . ,rt4 satisfying the closure condition 



jinH \-j 4 n 4 = 0. (12) 

The function $7(ni, . . ,714) provides normalization to one of the intertwiner. The function cj> m i"' m 4 
is invariant under rotations of the four vectors n\ , . . , n^. In the following we always assume that 
this invariance has been fixed with a given choice of orientation 5 . 

Nodes of the spin network can be labelled with coherent intertwiners. In fact such states provide 
an overcomplete basis of Tlr 5 - Calling ■y j mi '" m4 the standard recoupling basis for intertwiners, we 
can define the coefficients 

$i(ni,..,n 4 ) =v? 1 - m4 $ mi ...mt(n 1 ,..,n 4 ) . (13) 

We define a coherent spin network j a i,,& a ) as the state labelled by ten spins jab and 4x5 
normals n a b and given by the superposition 



|r 5 , jab, $a(n)) = ^ [Yl^ia(nab))\T 5 Jab,ia) • (14) 

a=l 

The expectation value of the metric operator on a coherent spin network is simply 

(r 5 , jab, <S>a\Ec- E c\ T 5Jab, $a) - ^jcajcb n ca ■ H cb (15) 

in the large spin limit. As a result we can choose the normals n a b so that the coherent spin network 
state is peaked on a given intrinsic geometry of S. 

Normals in different tetrahedra cannot be chosen independently if we want to peak on a Regge 
geometry [44] . The relation between normals is provided by the requirement that they are computed 
from the lengths of the edges of the triangulation A5. In fact, a state with generic normals 
(satisfying the closure condition (|12p ) is peaked on a discontinuous geometry. This fact can be seen 
in the following way: let us consider an edge of the triangulation A5; this edge is shared by three 
tetrahedra; for each tetrahedron we can compute the expectation value of the length operator for 
an edge in its boundary [42] : however in general the expectation value of the length of an edge seen 
from different tetrahedra will not be the same; this fact shows that the geometry is discontinuous. 
The requirement that the semiclassical state is peaked on a Regge geometry amounts to a number 
of relations between the labels n a b- In the case of the boundary of a Euclidean 4-simplex (excluding 
the 'rectangular' cases discussed in [M]), the normals turn out to be completely fixed once we give 
the areas of the ten triangles or equivalently the ten spins j a b, 

n ab = flabUcd) ■ (16) 

This assignment of normals guarantees that the geometry we are peaking on is Regge-like. In 
particular, in this paper we are interested in the case of a 4-simplex which is approximately regular. 



5 For instance we can fix this redundancy assuming that the sum j\fti + is in the positive z direction while 
the vector n\ x H2 in the positive y direction. Once chosen this orientation, the four unit- vectors Hx,..,fU (which 
satisfying the closure condition) depend only on two parameters. These two parameters can be chosen to be the 
dihedral angle cos #12 = Hi ■ ni and the twisting angle tan0 { i2)( 34 ) = %\ X X ^ \alxaf\ ■ 
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Sjab 

■7° 

perturbative expression for the normals solving the continuity condition is available: 



In this case the spins labelling the links are of the form j a b = jo + 5j a b with -4^ <C 1 and a 



n a b 



Uo + V) = n ab (h) + E » (a6)(a °#«i • ( 17 ) 



The coefficients computed in terms of the derivative of the normals n a b (for a given 

choice of orientation, see footnote [5]) with respect to the ten edge lengths, using the Jacobian of 
the transformation from the ten areas to the ten edge lengths of the 4-simplex. 

In the following we are interested in superpositions over spins of coherent spin networks. As 
coherent intertwiners are defined only up to a spin-dependent arbitrary phase, a choice is in order. 
We make the canonical choice of phases described in [35] . We briefly recall it here. Consider a 
non-degenerate Euclidean 4-simplex; two tetrahedra t a and tb are glued at the triangle f a b = fba- 
Now, two congruent triangles f a b and fb a m ^ 3 can be made to coincide via a unique rotation 
R a b £ SO (3) which, together with a translation, takes one outward-pointing normal to minus the 
other one, 

Rabflab = -n ba ■ (18) 

The canonical choice of phase for the spin coherent states \j a b,n a b) and Ijab^ba) entering the 
coherent intertwiners <3? a and is given by lifting the rotation R a b to a SU(2) transformation g a b 
and requiring that 

\j a b,n ba ) = D {j - b) (g ab ) J\j ab ,n a b) (19) 
where J : TCj — > TCj is the standard antilinear map for SU(2) representations defined by 

(e|(|a)<8 J\0)) = (/3|a) with \a),\0) e Hj (20) 

and (e| is the unique intertwiner in TLj ® TCj. In the following we will always work with coherent 
spin networks |Ts , j a b, <& a (n(j))) satisfying the continuity condition, and with the canonical choice 
for the arbitrary phases of coherent states. From now on we use the shorter notation \j, 3?(n)). 

Coherent spin networks are eigenstates of the diagonal components of the metric operator, 
namely the area operator for the triangles of A5. The extrinsic curvature to the manifold E 
measures the amount of change of the 4-normal to S, parallel transporting it along S. In a piecewise- 
flat context, the extrinsic curvature has support on triangles, that is it is zero everywhere except 
that on triangles. For a triangle f a b, the extrinsic curvature K a b is given by the angle between the 
4-normals Na and ./V^ to two tetrahedra t a and tb sharing the face f a b- As the extrinsic curvature 
is the momentum conjugate to the intrinsic geometry, we have that a semiclassical state cannot 
be an eigenstate of the area as it would not be peaked on a given extrinsic curvature. In order 
to define a state peaked both on intrinsic and extrinsic geometry, we consider a superposition of 
coherent spin networks, 

l*o> = E ^0 UUHft)), (21) 

j a b 

with coefficients ipj 6 (j) given by a gaussian times a phase, 

^0 (i) = ^ exp ( - £ a^M M_Jocd\ exp f z £ rQ b {jab _ jQab) \ (22) 
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As we are interested in a boundary configuration peaked on the geometry of a regular 4-simplex, 
we choose all the background spins to be equal, joaft = Jo- Later we will consider an asymptotic 
expansion for large jo- The phases 4>Q b are also chosen to be equal. The extrinsic curvature at the 
face f a b in a regular 4-simplex is K ao = arccosiV a • N = arccos(— |). In Ashtekar-Barbero variables 
(Eq, A)) we have 

00 = <f>o = lK ab = 7 arccos(-l/4) . (23) 

The 10 x 10 matrix 

a {ab)(cd) ig 

assumed to be complex with positive definite real part. Moreover 
we require that it has the symmetries of a regular 4-simplex. We introduce the matrices p( ab ^ cd ^ 
with k = 0, 1, 2 defined as 

p(ab)(cd) _ ^ .£ _ zero otherwise, (24) 

or a permutation of it and zero otherwise, (25) 
p(ab)(cd) _ ^ _^ otherwise. (26) 

Their meaning is simple: a couple (ab) identifies a link of the graph T§; two links can be either 
coincident, or touching at a node, or disjoint. The matrices pjf" b ^ cd ^ correspond to these three 
different cases. Using the basis p^ ab ^ cd ^ we can write the matrix a^ ab ^ cd ^ as 



a {ab){cd) = ST ni p(ab){cd) 



£ a k Pi ab)(cd) . (27) 



k=0 



As a result our ansatz for a semiclassical boundary state \^o) is labelled by a (large) half-integer 
jo and has only three complex free parameters, the numbers a^. 



3 The new spin foam dynamics 

The dynamics is implemented in terms of a spin foam functional {W\. Here we are interested 
in its components on the Hilbert space spanned by spin networks with graph T§. The sum over 
two-complexes can be implemented in terms of a formal perturbative expansion in the parameter 
A of a Group Field Theory [36]: 

(W\r 5 ,j ab ,i a ) = Y J ^ Na W(a) (28) 

The sum is over spinfoams (colored 2-complexes) whose boundary is the spin network (T^, j ao ,i a ), 
W(a) is the spinfoam amplitude 

w{a) n w f n w - ( 29 ) 

/Ccr vCcr 

where W v and Wt are the vertex and face amplitude respectively. The quantity N a in (|28p is 
the number of vertices in the spin foam a, therefore the formal expansion in A is in fact a vertex 
expansion. 
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The spin foam models we consider here are the EPRL 7 [19] and FK 7 [20] models. We restrict 
attention to < 7 < 1; in this case the two models coincide. The vertex amplitude is given by 

W v (j ab ,i a ) = {15j}(j a V+) imUab^a) Ufa-Uab) (30) 

a 

fc a •'a 

where the unbalanced spins j + ,j~ are 

j'ab = 7 ± ja&, 7 ± = — • ( 31 ) 

This relation puts restrictions 6 on the value of 7 and of j ab . The fusion coefficients f l +.+ (jab) are 
defined in [19] (see also [37]) and built out of the intertwiner v™ 1 '" mi in Hj 1 <g> • • ®Ttj 4 and the 
intertwiners v f . 4 in 7i -± (g> • • • ® H„±. Defining a map Y : H, — > H,+ (8> H,- with matrix 

'± Jl J4 J J J 

elements Y™ +m _ = (j + ,m + ; j~ ,m~\ Y\j, m) given by Clebsh-Gordan coefficients, we have that the 
fusion coefficients /*+ p+ are given by 

fia _ V ...V „"11'""14 „, m l--- m t ,. m l ••• m 4 fooN 

J i+i+ ~ I mim+m^ 1 m i m+ m - V i V i+ V i- " 

Indices are raised and lowered with the Wigner metric. 

Throughout this paper we will restrict attention to the lowest order in the vertex expansion. 
To this order, the boundary amplitude of a spin network state with graph T$ is given by 

(W\T 5 ,j ab ,i a ) = n(jab)W v (j ab ,i a ) , (33) 

i.e. it involves a single spin foam vertex. 

The function fj, is defined as fi(j) = X\. ab Wf ab (j). A natural choice for the face amplitude is 
Wf(j + ,D = (2j+ + l)(2j- + 1) = (1 - 7 2 )j 2 + 2j + 1. Other choices can be considered. We 
assume that n(\j a b) scales as A p for some p for large A. We will show in the following that, at the 
leading order in large jo, the LQG propagator ([3]) is in fact independent from the choice of face 
amplitude, namely from the function [J,(j). 



4 LQG propagator: integral formula 

In this section we define the LQG propagator and then provide an integral formula for it. The dy- 
namical expectation value of an operator O on the state \^q) is defined via the following expression 

= (jW . (34) 



6 Formula (|30[) is well-defined only for j half-integer. As a result, for a fixed value of 7, there are restrictions on 
the boundary spin j. For instance, if we choose 7 = 1/n with n integer, then we have that j has to be integer and 
j > n, i.e. j G {n, n + 1, n + 2, • • • }. 
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The geometric operator we are interested in is the metric operator E® ■ E\ discussed in section [2J 
We focus on the connected two-point correlation function on a semiclassical boundary state 

|^o)- It is defined as 

/~<abcd / T?a ipb jpc jpd \ / jpa j-ift \ i t?c jpd \ i'qcN 

We are interested in computing this quantity using the boundary state |^o) introduced in section 
[2] and the spin foam dynamics (|30f) . This is what we call the LQG propagator. As the boundary 
state is a superposition of coherent spin networks, the LQG propagator involves terms of the form 
(W\0\j,$(n)). Its explicit formula is 

r abcd _ gj Hj)(W\E--E b n E- m -Et\jMn)) _ S 7 j>(j)(W\E--E b n \jMn)) E.j m{W\E^Ej\jMn)) . , 

In the following two subsections we recall the integral formula for the amplitude of a coherent spin 
network (W\j, <£(«)) [32J, [34] , [35] and derive analogous integral expressions for the amplitude with 
metric operator insertions (W\E®-E b \j, <3?(n)) and (W\E%-E^ E^-E^lj, $(n)). 



4.1 Integral formula for the amplitude of a coherent spin network 

The boundary amplitude of a coherent spin network \j, <£(n)) admits an integral representation 
[32} I34j.[35]. Here we go through its derivation as we will use a similar technique in next section. 

The boundary amplitude (W\j, <3?(n)) can be written as an integral over five copies of SU(2) x 
SU(2) (with respect to the Haar measure): 

r 5 

{W\j ab ,^ a (n)) = ^[Yl^ a (n)){W\j ab ,ia) = Kj) J ]J d 9ad9a llP ab (g + ,g~). (37) 

i a a 1=1 ab 

The function P ab {g + , g~) is given by 

P a \g + ,g-) = (jab,-n ba \Y^D^ b )((g+)^g+)^ D ^((g-)^ . ( 38 ) 

where the map Y is defined in section [3l Using the factorization property of spin coherent states, 

Y\j,n) = |j + ,n) ® \r,n) , (39) 
we have that the function P ab (g + , g~) factorizes as 

P a \g + ,g-) = P ah+ {g + )P ab -(g-) (40) 

with 

P ab± = (jab^fibalD^iigtr'g^ljab^ab) = ((\, ~Uba | (<7a ± )" 1 9? I ~, "ofc) ) ^ ■ (41) 

In the last equality we have used (again) the factorization property of spin coherent states to 
exponentiate the spin j^ b . In the following we will drop the 1/2 in \^,n a b) an d write always \n a b) 
for the coherent state in the fundamental representation. 
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The final expression we get is 



r 5 

(W\j, $(n)> = Mi) / II d 9a^9a e S (42) 

J a=l 

where the "action" S is given by the sum S = S + + S~ , with 

S ± = Y, 2 ^g{-n ab \(gtr 1 9i\n b a) • (43) 



ab 



4.2 LQG operators as group integral insertions 

In this section we use a similar technique to derive integral expressions for the expectation value of 
metric operators. In particular we show that 

r 5 

(W\E%-E b n \j ab ,* a (n))=»(j) / Udg+dg- qt(g + ,g-)e s (44) 

J a=l 



and that 



5 

(W\E^.E b n E c m .Ei\j ab ^ a (n))=^j) / Hdg+dg- qf(g + , g~) q%(g+ , g~) e S (45) 

J a=l 

where we assume 7 n^m and a, b,c,d / n, m. The expression for the insertions q^{g + ,g^) in the 
integral is derived below. 

We start focusing on (E^-E^) in the case a / b. The metric field {E h a )i acts on a state \ j ab , m ab ) 
as 7 times the generator J« of SU(2). As a result we can introduce a quantity Qf 3 defined as 

Qt (g + ,g~) = {3a b ,-n b a\Y^D^\{gtT 1 gt)®D^\{g~r^ , (46) 

so that 

f 5 

(W\E na -E nb \j, $(n)> = J ds+ds" ^QrQf n /pC "(5 + ' ^) • (47) 



a=l cd 



The product n' i s over couples (cd) different from (na), (rife). Thanks to the invariance properties 
of the map Y, we have that 



Yjf\j ab , m ab ) = ( J? + + J? ~)Y\j ab , m ab ) . (48) 

Thus Qf b can be written as 

Qf> = Qf+ P ab - + P ab+ Qf~ (49) 

with 

Qi b± = 7 (ja^-^balD^dg^g^jf^j^Hab) . (50) 



7 Similar formulae can be found also in the remaining cases but are not needed for the calculation of the LQG 
propagator. 
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Now we show that Q" fe± is given by a function A^ b± linear in the spin j times the quantity p a6± 
defined in flU}, 

gab± = A ab± pab± _ ^-g 

The function A ab ^ is determined as follows. The generator J°" bziz of SU(2) in representation can 
be obtained as the derivative 



da 1 v v 11 



J, 



ab± 



(52) 



a'=0 



where the group element h{a) is defined via the canonical parametrization h{a) = exp(— ia l ^f) 
Therefore, we can write Q a as 





ct ! =0 



i>l (-™ba\{gt) l 9b h ^)\n a b)) 



a*=0 
±\-l„±| 



7 Jab (-^ba\(g a ) % ° l \nab){(-n b a\(g a ) 9 b \n a b)) 



Comparing expression (|53p with (|5ip and (|4ip . we find that ^4"" a± is given by 



/Jna 



("4n I (5a) ^nKna) 

A vectorial expression for A" a± can be given, introducing the rotation = D^(g^) 

-1 Rfi n*na Rq, ^an ^{Rn ^na X Rg, ^an) 



1 - (R a n an ) ■ (R n n r , 



Thanks to (|49p and (|5ip . we have that the expression for Qf 3 simplifies to 



with 



A 



A^ b p a ^ 



ab ^afe+ _|_ j^ab 



As a result, equation (|4T|) reduces to 

5 

(W|^-^|i,$(n))= / n d 5a + d5a ^^fII pCd ^ + '^ 



a=l cd 
na /info 



(53) 
(54) 

(55) 

(56) 
(57) 

(58) 



which is of the form (J42J) with the insertion A na ■ A . Therefore, comparing with equation (|44l) . 
we have that 

q a n b (g + ,g-) = A na -A nb (59) 

for a b. The case with a = b can be computed using a similar technique but the result is rather 
simple and expected, thus we just state it 



C(<? + ,<T) =7 2 jna(jna + l) 



(60) 
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As far as (E^-E^ E^-E^} a similar result can be found. In particular, for n / m and a, b,c,d / 
n,m the result is stated at the beginning of this section, equation (|45p . with the same expression 
for the insertion q c ^(g + , g~) as in equation ([59]) and equation (|60j) . 

Substituting (|44p - (|45p in ()36[) we obtain a new expression for the propagator in terms of group 
integrals: 

This expression with metric operators written as insertions in an integral is the starting point for 
the large jo asymptotic analysis of next section. 



5 LQG propagator: stationary phase approximation 

The correlation function (jBTj) depends on the scale jo fixed by the boundary state. We are interested 
in computing its asymptotic expansion for large jo- The technique we use is an (extended) stationary 
phase approximation of a multiple integral over both spins and group elements. In 15.11 we put 
expression (|61|) in a form to which this approximation can be applied. Then in 15.21 we recall a 
standard result in asymptotic analysis regarding connected two-point functions and in I5.3ti5.4l we 
apply it to our problem. 



5.1 The total action and the extended integral 

We introduce the "total action" defined as Stot = 1°8 ^ + & or m ore explicitly as 

Of- + — \ 1 (ab)(cd) Jab — JOab jcd — jocd ■ ST~^ , a b l ■ ■ \ 

Stot{jab,gZ,g a ) = ~os> a 7^= 7^= t Z^ ( P0 [Jab-JOab) 

+ S + (jab,gt) + S-(j ab ,g-) ■ (62) 

Notice that the action S + + S~ is a homogeneous function of the spins j a b therefore, rescaling the 
spins jo a b and j a f, by an interger A so that joab ~~ > ^JOab an d j a b — * ^Jab, we have that the total 
action goes to .Stot - > AStot- We recall also that — ► X 2 qf l b . In the large A limit, the sums over 
spins in expression (I6ip can be approximated with integrals over continuous spin variables 8 : 

H j d 5 g ± q ab e A5tot = J d w j d 5 g ± » qf e XStot + 0(X~ N ) VN>0. (63) 

3 

Moreover, notice that the action, the measure and the insertions in (|6ip are invariant under a SO (4) 
symmetry that makes an integration dg + dg~ redundant. We can factor out one 50(4) volume, 
e.g. putting gi = gi = 1, so that we end up with an integral over d 4 g ± = Y\ a= 2^9t ^9a ■ 



8 The remainder, i.e. the difference between the sum and the integral, can be estimated via Euler-Maclaurin 
summation formula. This approximation does not affect any finite order in the computation of the LQG propagator. 
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As a result we can re-write expression (161 1 in the following integral form 

abcd 4 |f]d^«e^ _ J d 10 j d V// gg e A ^ J d 10 j d V/x gff e A ^°* 

nm 1 /d 10 jdV/^ AS to t Jd 10 jdV/^e A ^ /d 10 jdV/^e A5 -t ^ lD4J 

To this expression we can apply the standard result stated in the following section. 

5.2 Asymptotic formula for connected two-point functions 

Consider the integral 

F(X) = [ dxf(x)e xs(x) (65) 



over a region of M. d , with S(x) and f(x) smooth complex-valued functions such that the real part 
of S is negative or vanishing, Re5 < 0. Assume also that the stationary points xq of S are isolated 
so that the Hessian at a stationary point H = S"{xq) is non-singular, detH ^ 0. Under these 
hypothesis an asymptotic expansion of the integral F for large A is available: it is an extension of 
the standard stationary phase approximation that takes into account the fact that the action S is 
complex [48 1 . A key role is played by critical points, i.e. stationary points xq for which the real 
part of the action vanishes, ReS^o) = 0. Here we assume that there is a unique critical point. 
Then the asymptotic expansion of F{\) for large A is given by 

/27r\ 2 pUndH \S(x ) / 11 \ 

F(A) = (t) VldetJfl ( /(X0) + xk^ H 1" + °) + °<*>) < 66) 

with /■•'■ = d 2 f / dx l dx : > and Ind-ff is the index 9 of the Hessian. The term D does not contain second 
derivatives of /, it contains only 10 f(xo) and fl(xo). Now we consider three smooth complex- valued 
functions g, h and /U. A connected 2-point function relative to the insertions g and h and w.r.t. the 
measure \i is defined as 

J dx fi(x) g(x)h(x) e xs ^ f dx /i(x) g(x) e xs ( x > f dx (j,(x) h(x) e xs ( x > 

/ dx n{x) e A5 ( x ) J dx fi(x) e xs ^ J dx/i(x)e xs ^ ' ^ ' 

Using (|66p it is straightforward to show that the (leading order) asymptotic formula for the con- 
nected 2-point function is simply 

G = ±(H-lyjg>(x )h' j (x ) . (69) 

Notice that both the measure function fi and the disconnected term D do not appear in the leading 
term of the connected 2-point function; nevertheless they are present in the higher orders (loop 



9 The index is denned in terms of the eigenvalues of of the Hessian as Indff = | J^ fc arg(/ifc) with — ^ < 
arg(ft fe ) < +f . 

10 More explicitly, the term D is given by 

D = f^R'JU^UH-'fiH- 1 )" 1 + ~f(x )I%j h (x )R^ nl (xo)(H-y™(H-y n (H- 1 ) hl (67) 

with R(x) = S(x) - S(x ) - \Hij(x )(x - x Y(x - x ) j ■ 
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contributions). The reason we are considering the quantity G, built from integrals of the type 
(|65p . is that the LQG propagator has exactly this form. Specifically, in sections 15.31 we determine 
the critical points of the total action, in 15.41 we compute the Hessian of the total action and the 
derivative of the insertions evaluated at the critical points, and in [7] we state our result. 



5.3 Critical points of the total action 

The real part of the total action is given by 

ReStot = - (Rea)( ab ^ 2ah ~Jl ah 3cd-JOcd + (7Q) 

ab,cd ^JOac VJOcd 

+ 2^ log 2 + zL-U lo g 2 • ( ^ 

Therefore, having assumed that the matrix a in the boundary state has positive definite real part, 
we have that the real part of the total action is negative or vanishing, ReStot — 0- ^ n particular 
the total action vanishes for the configuration of spins jab and group elements g^ satisfying 

jab = joab , (72) 
g± such that R±n ab (j) = -R b n ba (j) . (73) 

Now we study the stationary points of the total action and show that there is a unique stationary 
point for which ReStot vanishes. 

The analysis of stationary points of the action S + + S~ with respect to variations of the group 
variables g^ has been performed in full detail by Barrett et al. in [35]. Here we briefly summarize 
their result as they apply unchanged to the total action. We invite the reader to look at the original 
reference for a detailed derivation and a geometrical interpretation of the result. 

The requirement that the variation of the total action with respect to the group variables g^ 
vanishes, 5 9 S to t = 0, leads to the two sets of equations (respectively for the real and the imaginary 
part of the variation): 

\- .± Rgngb - R b n ha _ ^ , ± (Rtn ab ) X (Rfn ba ) _ 

^ a Jab l-(RUab)-( R b t n ba ) ' ^ a Jab l-(R^n ab )-(R^n ba ) " 1 ' 



When evaluated at the maximum point (|73h . these two sets of equations are trivially satisfied. 
Infact the normals n ab in the boundary state are chosen to satisfy the closure condition (|12p at 
each node. Therefore the critical points in the group variables are given by all the solutions of 
equation (|T3"1) . 

For normals n ab which define non-degenerate tetrahedra and satisfy the continuity condition 
(fTUj) . the equation R a n ab = —R b n ba admits two distinct sets of solutions, up to global rotations. 
These two sets are related by parity. The two sets can be lifted to 577(2). We call them g+ and 
g~ . Out of them, four classes of solutions for the couple {gt i9a) can be found. They are given by 

(at, 9a), (9a, gt), (gi,gt), {ga,ga)- (75) 
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The geometrical interpretation is the following. The couples (j a bn a b, jabn a b) are interpreted as the 
selfdual and anti-selfdual parts (with respect to some "time" direction, e.g. (0,0,0,1)) of area 
bivectors associated to triangles in 4-dimensions; since these bivectors are diagonal, they live in the 
3-dimensional subspace of M 4 orthogonal to the chosen "time" direction. Because of the closure 
condition (fT2l) . for a fixed n the four bivectors (j n an n a, jnatina) define an embedding of a tetrahedron 
in R 4 . The two group elements g£ and g~ of the action (|43p define an SO(4) element which rotates 
the "initial" tetrahedron. The system (173j) is a gluing condition between tetrahedra. The first 
two classes of solutions in (|75[) glue five tetrahedra into two Euclidean non-degenerate 4-simplices 
related by a reflection, while the second two classes correspond to degenerate configurations with 
the 4-simplex living in the three-dimensional plane orthogonal to the chosen "time" direction. 

The evaluation of the action S(j ao , g+ , g~) = S + (j ao , g£) + S~ (j a b, g~) on the four classes of 
critical points gives 

S(jab,gti9a) = +SR cgge (jab) , (76) 
S(jab,9ai9t) = ~ SReggetiab) , (77) 
S{jab,gi,g~i) = +7~ 1 S'Regge(ja£>) , (78) 
S(jab,9a^9a) = ~ l' 1 ^Regge (jab) , (79) 

where SR egge (j ao ) is Regge action for a single 4-simplex with triangle areas A ao = ^j a b and dihedral 
angles (fi a b(j) written in terms of the areas 

SReggeOafc) = ^ljab4>ab{j) ■ (80) 
ab 

Now we focus on stationarity of the total action with respect to variations of the spin labels j a b- 
We fix the group elements (gt -,9a) to belong to one of the four classes ([75]) . For the first class we 
find 

rggg^-^A, (si) 

(3i,9a) VJOabVJOcd 9j ab 

The quantity dSR egge /dj a b is 7 times the extrinsic curvature at the triangle f a b of the boundary of 
a 4-simplex with triangle areas A a b = 7 jab- As the phase <f>Q b in the boundary state is choosen to be 
exactly 7 times the extrinsic curvature, we have that equation (|8ip vanishes for j a b = joab- Notice 
that, besides being a stationary point, this is also a critical point of the total action as stated in 

On the other hand, if equation (|8ip is evaluated on group elements belonging to the classes 
(9a i9a)-> (9a i9a)i (9 a i9a)i we have that there is no cancellation of phases and therefore no sta- 
tionary point with respect to variations of spins. This is the feature of the phase of the boundary 
state: it selects a classical contribution to the asymptotics of a spin foam model, a fact first noticed 
by Rovelli for the Barrett-Crane model in [7j. 

5.4 Hessian of the total action and derivatives of the insertions 

Here we compute the Hessian matrix of the total action Stot sit the critical point 2 ah — joabi 
i9a->9a) = {9a 1 9a)' ^ e introduce a local chart of coordinates (p^,p~) in a neighborhood of the 



o as " 



dj 



ab 
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point (g£,g a ) on SU{2) x SU{2). The parametrization is denned as follows: we introduce 

gt(pt) = Hpt)9t (82) 



with h(Pa) = ^Jl — \p£ | 2 + ip^ • a. The vector p^ is assumed to be in a neighborhood of the origin, 
which corresponds to the critical point . We introduce also the notation 



n 



± 



(83) 



where is the rotation associated to the SU{2) group element g^. The bivectors {j a b n ab^ jab n ab) 
have the geometrical interpretation of area bivectors associated to the triangles of a 4-simplex with 
faces of area proportional to j a b- 

The Hessian matrix is obtained computing second derivatives of the total action with respect 
to jab, pt and p~ , and evaluating it at the point j a b = joab and p^ = 0. With this definitions we 
have that the (gauge-fixed) Hessian matrix is a (10 + 12 + 12) x (10 + 12 + 12) matrix (as it does 
not contain derivatives w.r.t. gf) and has the following structure: 



D tot — 



/ d 2 Stot 

djdj 



as 



d 2 S t , 



0l0xl2 Oioxl2 ^ 

d 2 S u 



0l2xl0 dp +dp+ °12xl2 
d 2 Stot 







12x10 







12x12 



For the non-vanishing entries we find 



d 2 s. 



tot 



p=0 



Q(ab){cd) 



tot 



if. 



± 

(ai)(bj) 



djabdjcd 

d 2 s< 



a 



djabdpi 



(ab)(cd) 



. 



(84) 



(85) 



p=0 



p=0 



'tot 



H, 



± 

(ai)(aj) 



d 2 Stot 



VJOabVJOcd + {S '^^ ' 

Kl ± 3*ab{^-nf b nl± + l e^n k ±) 



dp^dpt 



p=0 



p=0 



b+a 



% ab n ab ) ' 



(86) 
(87) 
(88) 



where we have defined the 10x10 matrix of second derivatives of the Regge action 

d 2 S Rc 



("SRcggc) (ab)(cd) 



icggc 



djabdjcd 



(89) 



]0ab 



16 



We report also the first derivatives of the insertion (g + ,g ) evaluated at the critical point: 

H 2 7 ± j0naj0nfe(^ - U na ■ U nb H± a + inta X «n&) ' 

-H 2 7 ± j0naj0nb{nta + "n&X 1 ~ ' ™«&) ' 

2 9(jna^na ' jnb n nb) 



dqf 




dpt 


p= 


dqf 






p= 


dqf 




djcd 


p= 



7 



dj 



cd 



(90) 
(91) 
(92) 



We recall that in all these expressions the normals n a b are functions of jab as explained in section [5J 
These expressions will be used in section [7] to compute the leading order of the LQG propagator. 



6 Expectation value of metric operators 

Before focusing on the LQG propagator, i.e. on the two-point function, here we briefly discuss 
the one-point function (E® ■ E%). Its meaning is the dynamical expectation value of the metric 
operator. The fact that it is non-vanishing provides the background for the propagator. Using the 
technique developed in the previous sections we can compute it at the leading order in the large 
spin expansion. We use the integral formula for the metric operator ()44|) -(H5l) and the stationary 
phase analysis of section [5] and find that the expectation value of the metric operator is simply 
given by the evaluation of the insertion q^(g + , g~) at the critical point 

(E-.E b n ) = qt(g+,g-) +O(j ). (93) 

For the diagonal components a = b we have that the insertion is simply given by = (jjna) 2 so 
that its evaluation at the critical point gives the area square of the triangle f na . For the off-diagonal 
components we have that q^ = ■ A n where is given in equation (|55p . Its evaluation at the 
critical point can be easily found using equation (|73|) in expression (i55j) . We find 



joab,9 + ,g~ 



+ A r ' 

jdab,g + 



30ab,g 



ijona^na (jo ) + ijona^na (jo ) = ijona^na (jo ) (94) 



so that A na at the critical point evaluates to the classical value £^ cl = 7Jona^na(jo)) the normal to 
the face a of the tetrahedron n (normalized to the area of the face) . It is the classical counterpart 
of the operator (E^) 1 . Therefore we have that at the leading order the expectation value of the 
off-diagonal components is given by the dihedral angle between two faces of a tetrahedron 

i E n ■ E n) = E nc\ ' E ncl + C(jo) 

= I 2 j0naj0nb n na {jo) • ^nb(jo) + O(j ) . (95) 

They have the expected geometrical meaning. We observe that the same quantities computed 
with the Barrett-Crane spinfoam dynamics do not show the right behavior when the off-diagonal 
components of the metric operator are considered. 
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Using the same technique we can evaluate the leading order of the two-point function. We have 
that 

( E n- E n E m' E m) = E ncV E nc\ E mcV E md + ^Oo) • ( 96 ) 

The quantity we are specifically interested in in this paper is the connected two-point function. 
It is of order 0{jq), therefore it requires the next-to- leading orders in equations (|95j) and (|96p . 
Such orders depend on the measure fJ-(j). However in the computation of the connected part, these 
contributions cancel. The technique we use in next section for the calculation of the connected 
two-point function is the one introduced in section (|5.2p and captures directly the leading order. 



7 LQG propagator: the leading order 

We have defined the LQG propagator as the connected two-point function G^ff = (E% ■ E h n E^ ■ 
E^ l ) — {E^-E^ l ) (E^-E^j). Using the integral formula ([4"I ]) -(|4"5" |) and the result (|69|) for the asymptotics 
of connected two-point functions, we can compute the LQG propagator in terms of (the inverse 
of) the Hessian of the total action and of the derivative of the metric operator insertions at the 
critical point. These two ingredients are computed in section 15.41 Using them, we find that the 
LQG propagator is given by 



dqf dqt , 

^nm W - 2^ V (P9)(«) /)»• fin + 



'(pq)(rs) P)A p>A 
p,q,r,s ° J P1 ° Jt 



rhi^ ){n){Sk) d P l + d P k s ++[H hrQWofrQrf- 

+ 0(jl) (97) 

where all the terms appearing in this expression are defined in section [531 From this expression we 
can extract the dependence on the boundary spin jo and on the Immirzi parameter 7. We notice 
that the combinations 

nabcd _ 1 n -l dtfn / Q8 \ 

7 J p< g, r <s °0p9 °Jrs 

1 V- V- ( 1 d 1n dq% , 1 dqf dqf 



yabcd _ 1 \^ ( _( rT+\-l Uq n U( im , 1 (tt-\-1 u( in u( im \ (QQ \ 

X nm ~ ^ [ 1+ ( H )(ri)(sk) g p k + + ^ W )(ri)(.fc) Q p i- g p k- ) > 

7", S — Z. I. tZ — _L 



dqf 


dqt 

dp k s + 


1 

+ — 


(H-) 


d P i + 


7 


dqf 




1 




dp l r + 


dp k s + 


1" 



Y ^ d = * y y fl-(F + )7 1 v ^ - — (H-)7\, 0<ti \ (10Q) 



are in fact independent from jo and from 7. In terms of these quantities we have that the LQG 
propagator has the following structure 

Gf^{a) = ( 7 jo) 3 (K b m d (a) + iX^f + 7 2 ^) + O(f ) (101) 
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abed j^abed 



where the dependence on 7 and on jo has been made explicit now. The matrices R nrr , . 
and Y^f 1 can be evaluated algebraically. Only the matrix R^f depends on the three parameters 
ao, ai, 012 appearing on the boundary state. We find that 



R 



abed 




\ 



where 



ci = 4/3i , c 2 = 4/? 2 



C3 



(2/3 - m + 3/3 2 ) 



C4 
and 11 

A) 
ft 



1 



(8/3 - 12ft) , c 5 = - (49/3 - 93/3i + 48/3 2 ) , c 6 



(23ft - 42 A + 15ft) 



1 

To 
1 

30 

1 

30 



1 



a + 6ai + 3a 2 + -8a - 8a% + 16a 2 + iVl5 a - 2a\ + a 2 + iVTE 

3 16 5 
+ 1= + 



32 



a + 6«i + 3a 2 -8a - 8ai + 16a 2 + iy/15 a - 2a% + a 2 + i\J\L 
3 64 5 



a + 6ai + 3a 2 -8a - 8«i + 16a 2 + a - 2ai + a 2 + £a/15 

The matrices X^ff and Y^ d turn out to be proportional 



X 



abed 



7 

36' 



zabed 



Y' 



abed 



36 



yabed 



(102) 



(103) 
(104) 

(107) 
(108) 
(109) 



(110) 



We notice that the inverse of the matrix Q( a b)(cd) can be written in terms of the parameters /3k using the formalism 
27|) introduced in section [2l 

2 

k = 



(Q ) 



(105) 



The matrix Q( a b)(cd) is denned in equation (|86|l and is given by 



Q 



(ab)(cd) 



70 



(ai>)(cd) 



(106) 



with ft = —7 a/| , hi 
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with the matrix Z^ d given by 



7abcd 







-1 


e 3 


-1 





e-*i 




e s 








e 3 


-1 


e 3 







-1 


e 3 









-1 




-1 







e s e 


—i— 

1 3 












(: 



















e*3 




e 3 







e 3 


-1 





e 3 
e ! 3 







— i- 
e 3 


-1 


— »'ZE 

e 3 





i- 

e 3 


-1 
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e 3 








i- 

e 3 


— i- 

e 3 


i- 

e 3 





-1 


e *3 


-1 







This is the main result of the paper, the scaling and the tensorial structure of metric correlations 
in LQG. In the following we collect some remarks on this result: 

- The LQG propagator scales as j'q, as expected for correlations of objects with dimensions of 
area square, E% ■ E h n ~ (7Jo) 2 - 

- The off-diagonal components are not suppressed as happened for the Barrett-Crane model 
[121 n~3] and have the same scaling as the diagonal ones. 



The contribution R^m i n (|10ip matches exactly with the matrix of correlations of areas and 
angles computed in perturbative quantum Regge calculus with a boundary state as done in 

D3H- 

On the other hand, the '7-terms' in (jlOip . jX^f - + j 2 Y^ d , are new and proper of the spin 
foam model. They come from SU{2) x SU(2) "group" fluctuations. They don't contribute to 
correlations, nor to area-angle correlations. On the other hand, their contribution 
to angle-angle correlations is non-trivial. 

In the limit 7—^0 and jo — > 00 with 7J0 = const = Aq, only the Regge contribution survives. 
It is interesting to notice that the same limit was considered in [39] in the context of loop 
quantum cosmology. 



The '7-terms' have an interesting feature that we now describe. Let us focus on the tensorial 
components G4 = G^ 4 ^ 45 ^ and G5 = G^ . They are related by a permutation of the 
vertices 4 and 5, keeping the other three vertices fixed. The 'Regge-term' is invariant under 
this permutation, i?^ 45 = ^12 . On the other hand the '7-terms' are not. In particular 
we have that 



7X (35)(45) + 7 2 Fi (35)(45) = ^ ^ X MW + ^JWj . (n2) 

It would be interesting to identify the origin of the phase We notice that the permutation 
of the vertex 4 with the vertex 5 of the boundary spin network corresponds to a parity 
transformation of the four-simplex. In this sense the '7-terms' are parity violating. 

In next section we investigate the relation of the result found with the graviton propagator 
computed in perturbative quantum field theory. 
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8 Comparison with perturbative quantum gravity 



The motivation for studying the LQG propagator comes from the fact that it probes a regime of 
the theory where predictions can be compared to the ones obtained perturbatively in a quantum 
field theory of gravitons on flat space [H [6] . Therefore it is interesting to investigate this relation 
already at the preliminary level of a single spin foam vertex studied in this paper. In this section 
we investigate this relation within the setting discussed in [12\ [131 E] ■ 

In perturbative quantum gravity 12 , the graviton propagator in the harmonic gauge is given by 



{huj,(x)hp(,(y)) 



-1 



2\x — y 



2 i^flp^VCr ^fjLCF^Up fi/wfi. 



(113) 



Correlations of geometrical quantities can be computed perturbatively in terms of the graviton 
propagator. For instance the angle at a point x n between two intersecting surfaces f na and f n b is 
given by 13 

q a n = g„v{x n )g pa {x n )BZ{x n )B v n l{x n ) (114) 

where B^ v is the bivector associated to the surface 14 . As a result the angle fluctuation can be 
written in terms of the graviton field 



ab\fiu 



(116) 



where we have defined the tensor (T^ 6 )^ 
is simply given by 



28 pa BZ{x n )B v n l(x n ). The angle correlation (G^%ft 



(G±%t = (v(^)v(^)) cc 6 r (?^r • (H7) 

In particular, this quantity can be computed for couples of surfaces identified by triangles of area 
Aq living on the boundary of a regular Euclidean 4-simplex. This quantity has been computed in 
|13j and we report it here for reference, 



G 



abed 




qft 18y/3 x 512 



Here we consider the Euclidean case. 
3 We thank E. Alesci for a discussion on this point. 



4 
6 



6 
6 

-16 )) 



14 To be more specific, we consider local coordinates (p , a 2 ) for a surface t and call t M (<r) its embedding in the Ad 
manifold. The bivector B^ u (x) is defined as 



dt" 



(115) 
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The question we want to answer here is if the quantity (G^Oqft and the leading order of the LQG 
propagator given by equation (jlOip can match. As we can identify jjo with the area Aq, we have 
that the two have the same scaling. The non-trivial part of the matching is the tensorial structure. 
Despite the fact that we have 9x9 tensorial components, only six of them are independent as the 
others are related by symmetries of the configuration we are considering. On the other hand the 
semiclassical boundary state |\&o) we used in the LQG calculation has only three free parameters, 
ao, «i, Q.2- Therefore we can ask if there is a choice of these 3 parameters such that we can satisfy 
the 6 independent equations given by the matching condition 

G»a)) iqg =((?tf) qft . (118) 

We find that a solution in terms of the parameters aj~ can be found only in the limit of vanishing 
Immirzi parameter, keeping constant the product jjo = Ao- In this limit we find a unique solution 
for otk given by 

1 



a 



a i 



|00 -(495616\/3 -45\/l5i) , (119) 
^(-299008^ + 35\/l5 i) , (120) 



a 2 = — (31744\/3 - bVl5i) . (121) 
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Therefore the matching condition (|118p can be satisfied, at least in the specific limit considered. 
Having found a non-trivial solution, it is interesting to study the real part of the matrix a^ ab ^ cd ^ 
in order to determine if it is positive definite. Its eigenvalues (with the associated degeneracy) are 

A 5 = 9216\/3 , deg = 5 , (122) 
d eg = 4, (123) 

Al= _12^, deg = 1 . {124) 

5 

We notice that all the eigenvalues are positive except one, Ai. The corresponding eigenvector 
represents conformal rescalings of the boundary state, jo a j, — ► Xjoab- It would be interesting to 
determine its origin and to understand how the result depends on the choice of gauge made for the 
graviton propagator (j!13j) . 



9 Conclusions 

In this paper we have studied correlation functions of metric operators in loop quantum gravity. 
The analysis presented involves two distinct ingredients: 

• The first is a setting for defining correlation functions. The setting is the boundary amplitude 
formalism. It involves a boundary semiclassical state \^o) which identifies the regime of 
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interest, loop quantum gravity operators E® ■ E n which probe the quantum geometry on the 
boundary, a spin foam model {W\ which implements the dynamics. The formalism allows to 
define semiclassical correlation functions in a background-independent context. 

• The second ingredient consists in an approximation scheme applied to the quantity defined 
above. It involves a vertex expansion and a large spin expansion. It allows to estimate the 
correlation functions explicitly. The explicit result can then be compared to the graviton 
propagator of perturbative quantum gravity. In this paper we focused on the lowest order in 
the vertex expansion and the leading order in the large spin expansion. 

The results found in the paper can be summarized as follows: 

1. In section [2] we have introduced a semiclassical state |^o) peaked on the intrinsic and the 
extrinsic geometry of the boundary of a regular Euclidean 4-simplex. The technique used to 
build this state is the following: (i) we use the coherent intertwiners introduced in [271 [34] 
to define coherent spin networks as in [35]; (ii) we choose the normals labelling intertwiners 
so that they are compatible with a simplicial 3-geometry (|16p . This addresses the issue 
of discontinuous lengths identified in [42] ; (iii) then we take a gaussian superposition over 
coherent spin networks in order to peak on extrinsic curvature as in [8] . This state is an 
improvement of the ansatz used in p21 [13j E] , as it depends only on the three free parameters 
ao, ai, ct2- 

2. In section U] we have defined expectation values of geometric observables on a semiclassical 
state. The LQG propagator is defined in equation (|35j) as a connected correlation function 
for the product of two metric operators 

/~iabcd / jpa jpb rpc jpd \ / rpa rpb\ / jpc jpd \ f-IOK\ 

This is the object that in principle can be compared to the graviton propagator on flat 
space: the background is coded in the expectation value of the geometric operators and the 
propagator measures correlations of fluctuations over this background. 

3. In section 14.21 we have introduced a technique which allows to write LQG metric operators 
as insertions in a 50(4) group integral. It can be interpreted as the covariant version of the 
LQG operators. The formalism works for arbitrary fixed triangulation. Having an integral 
formula for expectation values and correlations of metric operators allows to formulate the 
large spin expansion as a stationary phase approximation. The problem is studied in detail 
restricting attention to the lowest order in the vertex expansion, i.e. at the single-vertex level. 

4. The analysis of the large spin asymptotics is performed in section [5j The technique used is 
the one introduced by Barrett et al in [35]. There, the large spin asymptotics of the boundary 
amplitude of a coherent spin network is studied and four distinct critical points are found to 
contribute to the asymptotics. Two of them are related to different orientations of a 4-simplex. 
The other two come from selfdual configurations. Here, our boundary state is peaked also 
on extrinsic curvature. The feature of this boundary state is that it selects only one of the 
critical points, extracting expi5R e g ge from the asymptotics of the EPRL spin foam vertex. 
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This is a realization of the mechanism first identified by Rovelli in [7j for the Barrett-Crane 
model. 

6. In [6] we compute expectation values of LQG metric operators at leading order and find that 
they reproduce the intrinsic geometry of the boundary of a regular 4-simplex. 

7. Computing correlations of geometric operators requires going beyond the leading order in the 
large spin expansion. In section I5T21 we derive a formula for computing directly the connected 
two-point correlation function to the lowest non-trivial order in the large spin expansion. The 
formula is used in section 15.41 

8. The result of the calculation, the LQG propagator, is presented in section [7J We find that the 
result is the sum of two terms: a "Regge term" and a "7-term". The Regge term coincides 
with the correlations of areas and angles computed in Regge calculus with a boundary state 
|10j . It comes from correlations of fluctuations of the spin variables and depends on the 
parameters of the boundary state. The "7-term" comes from fluctuations of the £0(4) 
group variables. An explicit algebraic calculation of the tensorial components of the LQG 
propagator is presented. 

9. The LQG propagator can be compared to the graviton propagator. This is done in section [8l 
We find that the LQG propagator has the correct scaling behaviour. The three parameters 
cifc appearing in the semiclassical boundary state can be chosen so that the tensorial structure 
of the LQG propagator matches with the one of the graviton propagator. The matching is 
obtained in the limit 7^0 with 7^0 fixed. 

Now we would like to put these results in perspective with respect to the problem of extracting the 
low energy regime of loop quantum gravity and spin foams (see in particular [50J ) . 

Deriving the LQG propagator at the level of a single spin foam vertex is certainly only a first 
step. Within the setting of a vertex expansion, an analysis of the LQG propagator for a finite 
number of spinfoam vertices is needed. Some of the techniques developed in this paper generalize 
to this more general case. In particular superpositions of coherent spin networks can be used to build 
semiclassical states peaked on the intrinsic and the extrinsic curvature of an arbitrary boundary 
Regge geometry. Moreover, the expression of the LQG metric operator in terms of £0(4) group 
integrals presented in this paper works for an arbitrary number of spin foam vertices and allows to 
derive an integral representation of the LQG propagator in the general case, analogous to the one 
of [32] but with non-trivial insertions. This representation is the appropriate one for the analysis 
of the large spin asymptotics along the lines discussed for Regge calculus in |51j. The non-trivial 
question which needs to be answered then is if the semiclassical boundary state is able to enforce 
semiclassicality in the bulk. Another feature identified in this paper which appears to be general is 
that, besides the expected Regge contribution, correlations of LQG metric operators have a non- 
Regge contribution which is proper of the spin foam model. It would be interesting to investigate 
if this contribution propagates when more than a single spin foam vertex is considered. 
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